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Abstract 

In nuclear matter, the pion weak decay constant is separated into the two components ft, fs 
corresponding to the time and space components of the axial-vector current. Using QCD sum 
rules, we compute the two decay constants from the pseudoscalar-axial vector correlation function 
in the matter i f d'^x e^^''^' {p\T[d{x)ij^u{x) u{0)'y^^^d{0)]\p) . It is found that the sum rule for ft 



lO ' satisfies the in-medium Gell-Mann— Oakes— Renner (GOR) relation precisely while the fs sum rule 

does not. The fs sum rule contains the non- negligible contribution from the dimension 5 condensate 

CN ; 

^ ■ {qiD^iD^q) N + \ {qgs'^ • Qq)n in addition to the in-medium quark condensate. Using standard set 

, of QCD parameters and ignoring the in-medium change of the pion mass, we obtain ft = 105 

I MeV at the nuclear saturation density. The prediction for fg depends on values of the dimension 



5 condensate and on the Borel mass. However, the OPE constrains that fs/ ft ^ 1) which does not 
agree with the prediction from the in-medium chiral perturbation theory. Depending on the value 



Q^' of the dimension 5 condensate, fs at the saturation density is found to be in the range 112 ~ 134 

^ ■ MeV at the Borel mass ~ 1 GeV^. 

> : 
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I. INTRODUCTION 



The pion weak decay in nuclear matter is expected to be different from its decay in free 
space. In liadronic models, the in-medium shift of the decay constant occurs through the 
scalar meson exchange with background nucleons as well as through isobar-hole (or nucleon- 
hole) excitations that intermediates a pion and the axial- vector current The decay 
constant is directly connected to the renormalization of the induced pseudoscalar coupling, 
which therefore affects muon captures in nuclei ^, ^, ^ . 

One peculiar feature of the in-medium decay constant is its separation into the time and 
space components in the matter [0, ^. As we will see below, the separation seems to be 
a natural consequence when a nuclear matrix element of the axial-vector current involving 
a pion state is considered in defining pion decay constants. The two decay constants are 
scaled differently with the nuclear density and the induced pseudoscalar constant is separated 
correspondingly Interference between the two decay components may be important for 
explaining the muon capture data in nuclei. Therefore, it is necessary to determine the decay 
constants in the matter. Indeed, Kirchbach and Wirzba ||^ calculated the in-medium decay 
constants using the in- medium chiral perturbation theory developed in Ref.p. It is found 
that, while the time component is quenched moderately, the space component is quenched 
somewhat substantially. At the nuclear saturation density, its value is only 20 % of the free 
decay constant, which seems to be a too much quenching. 

Another important motivation for studying the in-medium decay constants is in con- 
nection with the so-called Brown- Rho scenario [Q. In this scenario, the in-medium decay 
constant is a scaling parameter that governs the in-medium reduction of hadron masses in a 
very simple manner. The quenching of the decay constant can be regarded as a signature of 
a partial restoration of chiral symmetry. To estimate the degree of the restoration, it is im- 
portant to calculate the in-medium decay constant model-independently. Such a calculation 
can be done using QCD sum rules 0, |^. 

In this work, we perform QCD sum rule calculation of the in-medium decay constant 
using the pseudoscalar-axial vector (PA) correlation function in nuclear matter. The in- 
medium chiral perturbation theory p[ will be used to represent the phenomenological side 
of the correlation function. In this approach, the low-lying intermediate state involves an in- 
medium quasi-pion whose pole is well separated from higher resonance poles. As we will see. 
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the sum rules constructed from this correlator satisfy in some limits the constraints coming 
purely from PCAC. Therefore, this correlation function may be relevant for studying pion 
properties in the matter. 

In the QCD side, we will use the finite density QCD sum rule techniques developed in 
Ref. 1^, |TU], |TT|. The finite density QCD sum rule has been used for investigating the self- 
energies of the nucleon embedded in the matter |r^, |IT] . Similar sum rule has been applied 
to the calculation of vector meson masses at finite density [l^ or temperature |T^. It can 



provide a useful constraint for various effective model predictions |]14|. A crucial element in 
the construction of these sum rules is the linear density approximation, which allows to write 
nuclear matrix elements in terms of nucleon matrix elements plus their vacuum expectation 
values. In-medium sum rules therefore are valid as long as the density is low. Following 
these in-medium frameworks, we construct sum rules for the PA correlation function. 

This paper is organized as follow. In Section II, we discuss a subtlety in the definition 
of the in-medium decay constants. The isospin symmetry must be assumed in the usual 



separation of the time and space components of the decay constant. In Section |T|, we 
construct the phenomenological side of the PA correlator by following the in-medium chiral 
perturbation theory. We then present a few constraints that have to be satisfied in some 
limits in Section |V]. The OPE calculation is given in Section |V] and the sum rule analysis 
for the decay constants is presented in Section |V|. 

II. SEPARATION INTO TWO IN-MEDIUM DECAY CONSTANTS 

There is a subtlety in the definition of the pion decay constants in nuclear matter. In 
this section, we address this issue by defining the in-medium pion decay constant in close 
analogies with the one in vacuum. To do this, we first recall that in vacuum, the pion decay 



constant fj, is defined by the relation |T5 



{0\dYl5u\7r+{k)) = if^k^ , {U = 131 MeV) . (1) 

A straightforward generalization to the nuclear matter case is to write the following nuclear 
matrix element into the form ^ 

(p| J7^75w|vr+(A;)p) = {ifth, ^/.k) , (2) 
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where \p) denotes the ground state of the nuclear matter. This is a usual definition that 
can be found in literature but one can show that taking this definition for the in-medium 
pion decay constants is equivalent to assuming an exact isospin symmetry in the isospin 
symmetric matter. The appearance of the two decay constants, ft and /g, is also related to 
the isospin symmetry. 

To illustrate this point in detail, we notice that the nuclear matrix element can contain 
a term proportional to the four- velocity vector (with = 1) of the nuclear matter as 
well as a term proportional to the pion momentum k^. The part proportional to can be 
further separated into two parts, one that vanishes in the soft-pion limit and the other that 
does not. Hence, we may write the nuclear matrix element into the form, 

(p|(i7''75u|7r+(A;)p) = iak^" + ibu^iu ■ k) + zcm^ , (3) 

with some in-medium (dimensionful) constants a, h and c. The third term that does not 
vanish in the soft-pion limit is proportional to the isospin breaking in the medium. This 
can be easily demonstrated by applying the soft-pion theorem directly to the nuclear matrix 
element, 

(p|J7^75M|7r+(A;)p)|fc,^^o ~ {p\d-i^'d - u7^m|p) ~ cu^" . (4) 

Thus, c in Eq. (^) is proportional to the isospin breaking in the medium and it is zero either 
when the isospin symmetry is exact or when the nuclear density is zero. Also under the 
isospin symmetry, there are only two in-medium constants involved for the matrix element 
Eq. @. In the nuclear rest frame = (1, 0), we recover Eq.(^ with the identifications ^ 

a + b = ft; a = fs . (5) 

When the isospin symmetry is broken in the matter, one may quickly notice that the time 
component ft is not well defined. In this case, the time component of Eq. is iak^ +ibk^ +ic 
and it is not clear how to deal with the c term in defining ft. Thus, ft defined in Eq. (H) is 
only for the isospin symmetric case. The ambiguity in the isospin breaking case is intrinsic 

^ One can expect that a and b in Eq. (||) also contain isospin breaking parts. As the isospin symmetry is 
assumed throughout this work, we will neglect such terms as well as the terms belonging to c. Thus, a 
and b in this equation are slightly different from the ones that appear in Eq. (||) by small isospin breaking 
parts. 



4 



in the matter and it may be important when the isospin breaking of the in-medium decays 
is studied. Somewhat larger breaking can be expected from the kaon channel as the strange 
vector current is involved. This breaking in fact makes and K~ decays different. In 
future, it will be interesting to study how the breaking effect is accommodated systematically 
in meson decays in the matter. In this work, we will neglect the isospin breaking and focus 
on the part that can be represented by Eq.(|^). The isospin breaking terms in the pion sum 
rules however are suppressed by the small quark masses and the density. 

III. PSEUDOSCALAR AND AXIAL VECTOR CORRELATION FUNCTION 

In our QCD sum rule study of the in-medium pion decay constants, we consider the 
pseudoscalar-axial vector (PA) correlation function in the nuclear medium. 



where \p) denotes the ground state of the nuclear matter. Because of the separation in 
Eq. (|^), there are two independent components for this in-medium correlator, time and 
space components. We will consider the following limit of the two components 



so that the resulting functions in the nuclear rest frame depend on one scalar variable p^. In 
this limit, the OPE calculation can be performed at — » —large unambiguously, which is 
then straightforwardly matched with the hadronic representation that will be constructed in 
this section. We will follow Ref. in making a phenomenological ansatz of this correlation 
function. 

Following the in-medium PCAC the pseudoscalar current couples to an in-medium 
quasi-pion. Hence, the low-lying state that intermediates the PA correlator is expected to be 
|7r(A;)p). Of course, other states containing vr' and tt" can also intermediate the correlator but 
we will treat them as the continuum. Introducing a coupling strength of the pseudoscalar 
current to the in-medium quasi-pion {p\di'~^^u\pT{) = G* and using Eqs. (|) (with c = 0), we 
obtain the phenomenological side of the PA correlator 0], 






(7) 



,^ apf^ + hu^iu ■ p) 



(8) 
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Here we have introduced the in- medium pion mass m*. As we will take the limits given in 
Eq. (|^, a possible p dependence in m* does not affect our sum rule. PCAC tells us that in 
vacuum, G* ^ = /^m^/2mq. In the matter, we should have G* = /tm*^/2mq so that it 
satisfies a set of relations given in Ref. [Q. Putting Eq. into Eq. (||) and evaluating it in 
the nuclear rest frame u'^ = (1,0), we obtain the phenomeno logical ansatz for the time and 
space components of the correlator n*^ and W (j = 1, 2, 3) in terms of ft and fg [defined in 
Eq. ©], 



IV. CONSTRAINTS 

In order to make a reliable prediction for the decay constants from the PA correlation 
function, we need to make sure that the intermediate state is dominated by the in-medium 
quasi-pion state. This quasi-pion dominance should be checked before one makes some 
claims on the in-medium properties from this PA correlation function. Here, we present two 
stringent constraints that have to be satisfied for the decay constants. These constraints 
will be compared with our sum rules later in order to justify our approach. 

One constraint is the well-known Gell-Mann— Oakes— Renner (GOR) relation 

-4mg(gg)o = m^/2 , (10) 

where the subscript denotes the vacuum expectation value. This has to be satisfied in the 
0{mq) order whenever we take the limit p — * 0. When p = 0, we have ft = fs = f-w and the 
decay constant is defined by Eq. (|I]). By taking the divergence of the axial- vector current 
and using the free equation of motion, we obtain 

2m,(0|Jz75M|vr+(A:)) = />2 . (H) 

Corrections to this will be of the order (9(m^) or higher. The GOR relation Eq.(0) can be 
easily obtained from this by applying the soft-pion theorem {k^ 0). 

The other constraint is the in-medium GOR relation which will constrain ft- Starting 
from the defining equation Eq. (^, we take the divergence of the axial- vector current again 
assuming that the quasi-pion is on its mass-shell with k'^ = m*^. It gives 

2m,{p\dzj,u\7i-^{k)p) = ftmf - (/, - ft)k' . (12) 



By taking k slightly off the mass shell, using the reduction formula to convert the pion state 
into a field, and following the similar steps as in deriving the soft-pion theorem |T^, we 
obtain 

-^-^{p\[dz^,u,Qt]\p)=ftm:\ (13) 
It 

where is the axial charge corresponding to the charged pion. Note, the soft-pion limit 
has eliminated the dependence on fg. Using the anticommutation relation for quarks, we 
easily evaluate the commutator to obtain the in-medium GOR relation 

- 4mq(gg)p = mff^ . (14) 

Again, this is an exact relation at the order 0{mq). The same relation has been obtained 
from in-medium chiral perturbation theory ||^, ^ |T^. Note that in deriving these two 
constraints Eqs.(p!0|) (0), we need to have a pion (or quasi-pion) state in the matrix elements. 
It means, the quasi-pion dominance in the PA correlation function can be justified once the 
resulting sum rules satisfy these constraints. 

V. THE OPERATOR PRODUCT EXPANSION (OPE) 

The OPE calculation for the pseudoscalar-axial vector (PA) correlation function Eq. (P) is 
performed in this section. For a simplicity, we will ignore possible contributions from instan- 



tons in this PA correlation function. For technical details, we mostly follow Refs. ||T0|, |19 
Before we proceed to the in-medium OPE calculation, it may worth mentioning on how the 
OPE looks like in the twist expansion. The leading term in the expansion, the scalar part, 
must take the same form as the vacuum sum rule. In the vacuum case, we know that 
the correlation function must reproduce the GOR relation at the order 0{mq). Note, the 
phenomenological side Eq. (^ contains quark mass in the denominator. To reproduce the 
GOR relation, only OPE term that does not vanish in the limit niq must be the quark 
condensate. All others in the OPE must be zero in the limit — > because the GOR 
relation is an exact relation at the order 0{mq). Modification to this picture in the medium 
is driven by the higher twist operators, which do not need to be scaled with quark mass but 
they must be proportional to the nuclear density. 

Another thing to mention is that the isospin symmetry in the medium is assumed in 
this work. It means that tt^ or vr^ decays are assumed to be the same. Note, once the 



isospin breaking is introduced, the time component ft is not well-defined. Anyway, under 
the isospin symmetry, the OPE must be symmetric under u ^ d as it should be. 

The OPE of the correlator Eq.(|) can be calculated by using the standard techniques in 
the fixed-point gauge [^. To include the finite quark mass effectively, we express the 
correlator in the momentum space. Introducing an in-medium quark propagator Sj{k) 
(j = u, d) in the momentum space, we write Eq. as 

/d'^k 
^Tr[z5,(fc)z75^5.(fc-p)7^5] . (15) 

The trace is taken over the Dirac and color indices. The perturbative part can be calculated 
by using the free propagator 

^Sr{k)=^^^^ U = u,d) (16) 

in Eq. (|l^). Using standard techniques of the dimensional regularization and the Feynman 
parameterization, we readily compute the perturbative part, 

n^ert(p') = -^'P" I mg\n[-u{l - u)p' + mj (17) 

with TUq = [nid + mu)/2. 

Nonperturbative contributions involving the quark condensate can be obtained by 
replacing one propagator by the disconnected nonlocal quark condensate, iS^^{k) — »• 
J d^x e*'^'^(g"(a;)gf (0))p while the other propagator remains to be the free quark propa- 
gator. The subscript p indicates the in-medium expectation value. A Tayler expansion 



around x,, = leads to 10, 21 



I d\ e^^-(g:(a:)gf(0)), = -^(2vr)^{5-/5A(A:)+7f5^(fc)} (18) 



where 



B\k) = {q^^q)^S^^)(k) + {q^'D,q),-^5^'\k) 

A few technical remarks are in order. Note that the dimension of the correlator is even. 
In the phenomenological side Eq. (P), one dimension is taken up by so the rest has an odd 
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dimension. It means, the constants a and h in Eq. (P) pick up nonzero contributions only 
from odd dimensional operators in the OPE. Even dimensional operators are not symmetric 
under u ^ d and constitute the isospin breaking parts of the OPE which we ignore in 
this work. Then, this observation yields interesting consequences. When the part A{k) in 
Eq. (|I^) is taken for the disconnected quark propagator, then only the |^ part of the free 
quark propagator Eq. (^) survives through the trace in Eq. (0). In this case, the term like 
{qD^j_q)p constitutes even dimensional contributions, thus cannot contribute to our isospin 
symmetric decomposition of Eq. (H). Similarly when the part in Eq. ([l8|) is taken for the 
disconnected quark propagator, then the irij part of the free quark propagator gives nonzero 
contributions. Odd dimensional operators in (for example {q'^^q) p) combine with irij 
to form even dimensional operators and they can not contribute to Eq. (Rf). In- medium 



condensates that contribute to our sum rules can be decomposed into [|I0 

4 
3 
4 
3 



{qDpD^q)p 
{q-fpD^q)p 



Mu-Dfq)^ 



u ■ Dq)^ 



-{qD^q)p{u^,u^- g^^) 



(20) 
(21) 



Another nonperturbative contributions come from the quark-gluon mixing operator. To 
form such an operator, we need to take the quark propagator with one gluon line attached p| 



iScd{k) 



9. 



2 ^''(P 



m 



2^3 



(^ + mg)7'^(^ + m,)7"(^ + m, 



(22) 



where G^^ is the gluon field tensor, the SU(3) generators in color space with the nor- 
malization Tr{t^t^) = /2, c and d color indices. Taking the gluon tensor and moving it 
into the disconnected quark propagator, one can form a quark-gluon mixed condensate in 
the medium |Tn 



+ {gsqicF ■ Qq)p [<Jpui + i{up'^u - Mj^T/^)]"^ 

- 4 [{q{u ■ Dfq) + img{q'^u ■ Dq)] [a^y + 2z(m^7^ - u^7^)^]°^ | . (23) 

Here note, Qpu = Gp^t"^. 

Combining all these ingredients, we readily compute nonperturbative parts that can con- 
tribute up to dimension 5. After a lengthy but straightforward calculation, we obtain 



nonpert 



— 



2 + 3"^9 



pf^ — Ap ■ uu^ 



2i 



m 



2^2 



{P^ 



m 



2^3 



+ 



32z 



— {q{u ■ D) q)p - —{qgsCT ■ Qq) 



p ■ uu^ {p ■ uYp^ 



m 



2\2 



(p2 



m 



2^3 



16 



P ■ UU^ — p^ 1 4 4 i-jn \ P ' 



(24) 

-gy \f ■■"qj 

Here, we have not written down even dimensional operators belonging to the isospin break- 
ing. 

At the dimension 5, one can expect another nonperturbative contributions (with odd 
dimensional condensate) containing f^q{^G'^)p (or gluonic twist operators times m^). To 
compute this, we use quark propagators with one and two gluon attached [§] in the calcu- 
lation of Eq. (|l^). Somewhat lengthy calculations show that the Wilson coefficient of this 
term is infrared divergent when the quark mass goes to zero. In the case of the vacuum 
sum rule, Ref. |^ proposed a simple prescription to extract the gluonic Wilson coefficient. 
In our in-medium sum rule, further complications arise due to gluonic twist operators and 
the simple prescription proposed in Ref. needs to be further developed for the twist 
operators. However, all these gluonic contributions must contain the quark mass (to make 
dimension 5 condensate) and they should be small. They may however be important in the 
kaon channel as the strange quark mass is much larger than m„ or m^, which will be studied 
in future. 

Since we are constructing sum rules for the two scalar functions given in Eq. (P), we can 
write down from Eqs. (|l^ (p^ the time and space components of H'^ = H^^^^ + H: 



nonpert 



m 



the nuclear rest frame = (1, 0) 



H = hm — 

"P"" P^O tpQ 



H 



ope 



4vr2 

lim — r 

p^o ipi 

3 



du mg\n[—u{l — u)pl + m^] + 



2{qq)p 8mq{qHDoq)p-2ml{qq)f 



pI 



+ 



(P 



m 



2^2 



47r2 
32 



du mg\n[—u{l — u)pl + m^] + ^^^^^^ 



+ - 



{qiDoiDoq)p + -{qgsCr ■ Qq)^ 



Po 
1 



8mg{qUDoq)^ 
3 (pI- 



2^2 



2ml{qq)p 

{pI - m; 



2^2 



{j>l - m2)2 



(26) 



In obtaining this, we have neglected the terms whose dimension in the numerator is 7. 
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VI. THE IN-MEDIUM PION DECAY CONSTANTS 



Having calculated the OPE for the PA correlation function, we can easily construct QCD 
sum rules for the in-medium pion decay constants ft and fs- Taking imaginary part of Eqs. 
(P), (p^), (pB) and matching them within the duality region with a Borel weight e~^^^^^ , 



Jo ^ 



{l = t,s) 



(27) 



we obtain sum rules for the time and space components separately, 
*2 rSo 



m. 



2m. 



-ft'e 



4m2 



2{qq)pe~ 



+ 



2m, 



-ftfsC 



QUlq 

Sniq 
4^ 



2ml 
M2 



(28) 



So 



4m2 



4m? 



+ 



3M2 
32 
3M2 



2m' 



9 0/^^\ ^-m-l/M'^ 



2{qq)pe- 



M2 



9 /™\ ^-ml/M^ 



(29) 



The continuum threshold is denoted by Sq. The dependence on this parameter is expected 
to be small as the perturbative part is proportional to the quark mass. 

Various nuclear matrix elements appearing in the OPE can be written in terms of cor- 
responding nucleon matrix elements (denoted by the subscript 'W below) in the linear 
density approximation. In this approximation, the dimension 4 and 5 condensates can be 
written |TD] 



{qUDoq)p = p {qHDoq)N + -^{qq) 



p ' 



(30) 



{qiDoiDoq)p + -{qgsCr ■ Qq)p 
o 



P 



m„ 



+ '-f{qq)p. (31) 



{qiDoiDoq)N + -{q9sO- ■ Gq)N 
o 

These expressions are slightly different from the ones in Ref. [l^ by the quark-mass depen- 
dent terms but whose appearance can be easily understood from the equation of motion. 
Ref. neglected these small quark mass terms. It is interesting to see that, with these 
quark mass terms, the Wilson coefficient of the condensate m^(gg)p is delicately canceled 
away. Anyway, in this separation, the nucleon matrix elements should be understood as the 
ones in the chiral limit m 



^ 0. According to Ref. [T^ 
{qHDoq)N ~ 0.18 GeV , 



(32) 
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{qiDotDoq)N + liqOsCr ■ gq)^ ~ 0.08 - 0.3 GeV^ . (33) 

o 

The value of the first equation comes from the twist-2 quark distribution function in a 
nucleon. The value of the second line, the dimension 5 term, is not well-known. This is 
related to the twist-3 nucleon distribution function whose value however can not be reliably 
calculated. The bag model estimate gives 0.08 GeV^ ||lOl while the other standard estimate 



yields 0.3 GeV^ [10|. In our later analysis, we will take these two values as a window of this 
dimension 5 condensate. 

But it is important to note that this dimension 5 term can not he negative. To prove this 
statement, first note that (g^'scr ■ Qq)N = '^{qD'^q)N- Using further that = Dq — D^, one 
can rewrite 

13 1 
{qiDoiDoq)N + -^{q9sCr ■ gq)^ = -{qiDoiDoq)N + -{q{i'Dfq)N ■ (34) 

Note, iDf) and iY) are hermitian operators and therefore their square must be positive 
definite. This means that this dimension 5 condensate, as it contains the positive definite 
operators, must have the same sign as the positive quantity {qq)N- 

Using Eqs.(pil|), (PT|), we rewrite Eqs.(P5|), ( |2^ ) as (transferring the quark mass from the 
LHS to the RHS) 

mf/^^e-*^/^^^ = ^ r dse-^'''^'J^^^-AmM)pe-<"'' 

27r^ J^rnl V S 



<3 

o2 /A,f2 



+ -^P {q'iD,q)^e-<l'^' , (35) 



2vr^ ^4^2 V s 



{qUDoq)r,e-'^/'^ 



3M2 



64mq 



{qiDoiDoq)N + \ {qgsCT ■ gq)N 



e-'^V'"' . (36) 



Our two sum rules mainly differ by the dimension 5 condensate {qiDi^iDi^q) iv + 1 {qgs^ ■ Qq) n- 
Other terms that make the two sum rules different are suppressed by higher orders in small 
quark mass. 

As we have announced, let us see if these sum rules satisfy the constraints discussed in 
Sec. II. First in the limit, p — 0, we have ft = fs = fn, = "^vr- Then, one can easily 
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see that the two sum rules become identical to the vacuum pion sum rule that appears in 
Ref. [m, 

/>^e-^/^^ = 1^ r ds e-^/^'Jl^^-4m,{qq)e-V^' . (37) 

J4m2 V S 

This give a mathematical check for our OPE expression. Furthermore, at the order 0{mq), 
one can see that the GOR relation is precisely satisfied (Note, the perturbative part is an 
order O^m"^).). 

The additional check can be made by observing that the ft sum rule satisfies the in- 
medium GOR relation. Namely, in the sum rule for ft Eq. (^), dominant piece of the 
OPE comes from 'mq{qq)p. The other OPE are down by higher orders in quark mass. That 
is, at the order 0{mq) [or at the order in (9(m*^)], this sum rule precisely reproduces the 
in-medium GOR relation [See also Ref. [Q, |, p!7|].], 

m*Jft^ ^ -Amq{qq)p . (38) 

Such a simple relation can not be obtained from the space component sum rule Eq. (|36| ) due 
to the non- negligible dimension 5 term {qiDQiD^q) ^ + \ {qgs(T ■ Qq)N- Based on these two 
checks, the quasi-pion dominance in the PA correlation function can be justified because the 
two constraints require a quasi-pion in the intermediate state. 

The ratio of Eqs. (|35|) (|36|) gives us an interesting relation for ft and fs- Namely neglecting 
the terms of order 0{m^^) and higher, we find the approximate formula, 

ll _ (g^-Po^-Pog)jv + Wdsfy ■ Qq)N , . 

ft"" 3M2 (gg), " ^^^^ 

When p ^ 0, the ratio becomes the unity as it should be. Note, {qq) p is negative and the 
dimension 5 condensate in the numerator is, as we have emphasized, positive. Therefore, we 
should have 

fs>ff (40) 

One may suspect that the OPE does not converge fast enough up to dimension 5 and 
large contributions from higher dimensional operators can alter this inequality. But non- 
negligible higher dimensional operators (even if they exist) must be related to twist operators 
that scale with the density. The scalar (nontwist) operators with higher dimensions must 
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be proportional to 0{m?^ or higher orders because the scalar part of the sum rule at the 
order 0{mq) must satisfy the GOR relation. Higher twist operators are suppressed either in 
the twist expansion and by the high correlator momentum. For operators having the same 
twist but with higher spin indices, we do not have a systematic way to analyze them but 
their contributions are usually smaller than the contribution from lower spin operators |T^, 
TBI in in-medium sum rules. Even if there are large corrections from higher dimensional 



operators (though it is unlikely), the Borel window from which the physical parameter is to 
be extracted will be chosen such a way that their contributions are small. In this window, 
their contributions will be down by higher orders of the Borel mass (M^, and so on) and 
their effects affect only the region around the lower boundary of the window. Therefore, it 
is unlikely that such higher dimensional terms reverse the inequality. 

On the other hand, the in-medium chiral perturbation theory gives the ratio in terms of 
the chiral perturbation parameters(Note, = 131 MeV in our notation.) [Q], 

il. = h = 1 h f41') 

ft 1 + 4(C2+C3)P 4(C2+C3)P ■ ^ 

fir f TT 

With the parameter values, C2 = 0.28 fm and C3 = —0.55 fm given in Ref. |^, the ratio is 
0.28 at the nuclear saturation density p = 0.17 fm~'^, substantially smaller than the unity. 
If one compares with the free decay constant, the value seems to be suppressed too much, 
only 20 % of f-^- Even if we take into account an error in fs coming from a large uncertainty 
in the parameter C3 0, our finding is still inconsistent with this result. To be comparable 
with this result, the higher dimensional terms that we have not calculated need to be an 
order of the leading term in the OPE, which certainly is not realistic. 

We now include the higher order rriq corrections in our analysis and calculate the decay 
constants at the nuclear saturation density. In the phenomenological side, we have three 
parameters to be determined ft, fs and m* but we have only two sum rules. The best fitting 
method with a trial function of the form ae~^^^^^ may not reliably determine the exponent 
b. Further input may be needed. For simplicity, we take the in-medium pion mass to be its 
free value m* = m^r. This assumption is supported by Ref. ^ where the mass change 

in the medium is calculated to be about 10 MeV. It should be noted that this assumption 
does not affect one of our results Eq. (|40|) as the pion mass is canceled in the ratio. 
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The in-medium quark condensate in the hnear density approximation is given by [ITn, |25 



im), = {qq)o + P 1^ • (42) 

The vacuum expectation value is denoted by the subscript "0" . In our analysis, we will use 
{qq)o = —(225 MeV)'^ that corresponds to the quark mass rrig ~ 7 MeV p^, |2^. The nucleon 
sigma term is cxn ~ 45 MeVp8|. The values for the other nucleon matrix elements are given 



in Eq. (pSf) . The continuum threshold 5*0 will be fixed to the vr' mass but the sensitivity of 
our result on this choice is very small. 

In Figure |I], the solid lines show the time component of the decay constant ft with respect 
to the Borel mass. The upper one is the vacuum case and the lower one is from Eq. ( pS] ) at 
the nuclear saturation density p = 0.17 fm^^. The dashed lines are obtained from the GOR 
relation, ft from the OPE is only 2 % larger than the one from the GOR relation. It is clear 
from this figure that there is no dependence on the Borel mass. At the saturation density, 
we obtain 

ft = 105 MeV , (43) 

20 % lower than its free space value. It is not so different from the value of the in-medium 
chiral perturbation theory ft = 101 MeV 0. 

Figure shows the Borel curve for the ratio fs/ ft- The solid curves are from Eqs.(p5D 
(P^) and the dashed curves are from our approximate formula Eq.(P^), which are almost 



indistinguishable from the full curves. Therefore, higher corrections are negligible. In 
getting the upper two curves, we use the dimension 5 condensate {qiDoiDoq) m + \ {qgscy ■ 
Qq)N = 0.3 GeV^ while the lower two curves use 0.08 GeV^. The Borel curves are not so 
stable with respect to the Borel mass. Therefore, they are not so useful in extracting the 
precise value of the ratio. Higher dimensional operators that we have not included in our 
calculation may affect the Borel curves at lower mass region and stabilize the curves. Their 
contributions however should become smaller as we go to higher Borel masses. Anyway, 
because of this Borel instability, we are not able to predict a precise value for the ratio. As 
a rough estimate, at = 1 GeV^, we find the ratio fs/ ft = 1-28 from the upper curve and 
1.07 from the lower curve. The variation within 0.8 GeV^ < < 1.2 GeV^ is estimated 
only 3 % - 10 %. Using our result of ft = 105 MeV, the two extremes give the range 

112 MeV < < 134 MeV . (44) 
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Thus, in this rough estimate, the space component of the pion decay constant either shghtly 
increases or decreases from its vacuum value depending upon the value of the dimension 5 
condensate. Of course a more detailed analysis after including higher dimensional operators 
may be need to determine the precise value of fg. 

VII. SUMMARY 

In this work, we have constructed QCD sum rules for the pion decay constant in nuclear 
matter. The pseudoscalar-axial vector correlation function in the matter has been used for 
this purpose. The sum rule for ft is found to satisfy the in-medium GOR relation, which 
indicates that the method we used is reasonable for investigating the pion properties in 
the matter. The OPE for the space component contains the non-negligible dimension 5 
contribution in addition to the quark condensate. We have established a firm constraint 
fs > ft and it does not agree with the result from the in-medium chiral perturbation theory. 
Neglecting the in-medium shift of the pion mass, we have obtained ft — 105 MeV. Using 
this value and from a rough estimate, the space component of the decay constant is found 
to be in the range 112 MeV < fs < 134 MeV. 
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FIG. 1: The Borel curves for ft, the sohd hnes from its sum rule and the dashed hnes from the 
GOR relation. The upper two curves are the vacuum case and the lower two curves are obtained 
from the in-medium sum rule at the nuclear saturation density p = 0.17 fm~^. 
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FIG. 2: The Borel curves for the ratio jsl ft at the nuclear saturation density p = 0.17 fm ^. 
The upper curves are obtained with the dimension 5 condensate 0.3 GeV^ and the lower curves are 
obtained with 0.08 GeV^. The solid curves are from our sum rules and the almost indistinguishable 
dashed curves are from the approximate formula Eq. (|39|). 
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